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We demonstrate an iterative scheme for coupled-cluster properties calculations without truncating 
the dressed properties operator. For validation, magnetic dipole hyperfine constants of alkaline Earth 
ions are calculated with relativistic coupled-cluster and role of electron correlation examined. Then, 
a detailed analysis of the higher order terms is carried out. Based on the results, we arrive at 
an optimal form of the dressed operator. Which we recommend for properties calculations with 
relativistic coupled-cluster theory. 
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I. INTRODUCTION 

Coupled-cluster theory, first developed in nuclear 
many body physics [H H] , is considered one of the best 
many body theory. In recent times, it has been used 
with great success in nuclear [3^ , atomic IH [5] , molecular 
^ and condensed matter [7J calculations. In atoms it 
is equivalent to incorporating electron correlation effects 
to all order. It has been used extensively in precision 
atomic properties and structure calculations. These in- 
clude atomic electric dipole moments [HH], parity non- 
conservation [3] , hyperfine structure constants [Sj [TD] and 
electromagnetic transition properties [TTJ [12] . 

Despite the remarkable developments and numerous 
calculations based on relativistic coupled-cluster theory. 
Hitherto, a systematic analysis of the properties calcu- 
lations with coupled-cluster wave functions is lacking. 
This issue arises from the fact that, the expression for 
properties with coupled-cluster wave functions is a non 
terminating series. In this paper we demonstrate an iter- 
ative scheme to calculate properties without truncation. 
Such a study is essential and timely as precision atomic 
calculations, in several instances, complement precision 
atomic experiments. These have direct bearing on, to 
mention a few, fundamental physics and new technology. 

To test and validate the scheme we employ open shell 
coupled-cluster theory |T31 HH [TS] and calculate the mag- 
netic dipole hyperfine constants of alkaline Earth ions 
25Mg+, '*3Ca+, 87Sr+ and ^^'^Ba+ . We have selected 
these ions as these are potential candidates for ongo- 
ing or proposed novel experiments. In addition, there 
is a large variation in the role of electron correlation 
among the ions and states. The ground state hyper- 
fine constant of Mg"*" is well studied with ion trapping 
techniques [TB]. The clock states of, the next ion in 
the group, ^'^Ca''' were recently employed for high-fidelity 
entanglement |17j . A crucial step in quantum informa- 
tion processing. Then, single trapped ^''Sr^ is a suitable 
frequency standard |18j . These are application oriented 
precision experiments. The other fascinating prospect is 
observation of parity nonconservation in a single ^"^^Ba^ 
|19j . In all of these endeavours, hyperfine interaction is 
involved. For this reason, several theoretical calculations 



have examined the role of electron correlations to the hy- 
perfine constants of these ions. These provide a wealth 
of data for comparative study. In addition to magnetic 
hyperfine constant, we also compute the excitation ener- 
gies of the low lying states. This is to verify the quality 
of the single particle wave function we use. 

The paper is divided into seven sections. In the next 
section, that is Section. II, we give a brief description of 
single valence coupled-cluster theory. Then, Section. Ill 
is a short writeup on hyperfine interaction and how it 
is calculated with relativistic coupled-cluster. Section. IV 
forms the core of the paper, where we explain our it- 
erative scheme to calculate properties with relativistic 
coupled-cluster to all order. The details of the numerical 
methods and schemes used in the present work are pro- 
vided in Section. V. And then we present our results in 
Section. VI. Finally, in Section. VII we make concluding 
remarks, which may serve as guideline for any properties 
calculations with relativistic coupled-cluster theory. In 
the paper, all the calculations and mathematical expres- 
sions are in atomic units {e — h~ rrie = 1). 



II. SINGLE VALENCE COUPLED-CLUSTER 
THEORY 

For completeness and easy reference of the working 
equations, we provide a condensed overview of the single 
valence coupled-cluster theory. Readers are referred to 
Ref. [in] for a detailed exposition of the theory. In the 
Fock space coupled-cluster theory of single valence sys- 
tems, the correlated wave function is calculated in two 
steps. First, the cluster operators of the core electrons 
or the closed-shell part T is evaluated from the reference 
state |<i>o)- Second, the cluster operators of the valence 
shells S is evaluated and the reference state is 

\'^.)^al\^o)- (1) 

The coupled-cluster wave function of the open shell sys- 
tem is 

=e^+^|a>,). (2) 
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For single valence system ^ 1 + S, the higher order 
terms in the exponential do not contribute. Then 

1^,) =6^(1 + 5)1$,). (3) 

For an N electron atom, the cluster operators are 



JV-l N 

T — y Ti, and S = 



1=1 



(4) 



Here the summation index of the T is up to the iV — 1 
core electrons, where as S is up to N to include the va- 
lence electron. However, single and double are the most 
dominant, in coupled-cluster singles and doubles (CCSD) 
approximation T — Ti +T2 and = 5i -|- ^2 . In the sec- 
ond quantized representation, for the closed-shell part 



=Y1 *a(4>^a 
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and T2 



1 
2! 



a,b,p,q 



(5) 



Similarly, for the valence shell 

Si = ^slalay, and 5*2 = ^ sHalalaaa^. (6) 

P Oj.p.q 

Here, t..'. and s..'. are the cluster amplitudes. The in- 
dexes abc . . . {pgr . . .) represent core (virtual) states and 
vwx . . . represent valence states. The operators Ti (Si ) 
and T2 {S2) give single and double replacements after op- 
erating on the closed- (open-) shell reference states. The 
diagrammatic representation of S are shown in Figjl] 
The atomic state |^^) satisfy the eigen value equation 



(7) 



where H is the atomic Hamiltonian and Ey is the exact 
eigen energy of the atomic state. Applying e^'^ on the 
above equation, we get 

H{l + S)\<i>y)=Ey{l + S)\<i>y), (8) 

where 

1 



=n 



H = H+{HT}+^{HTr}+^{HTTT}+^{HTTTT}, 

(9) 

is the dressed Hamiltonian and {• • • } denotes normal or- 
dering of the operators and {A ■ ■ ■ B} represents contrac- 
tion between two operators A and B. The cluster ampli- 
tude equations of the singles and doubles are obtained af- 
ter projecting Eq.(|8]) on singly and doubly replaced states 
(<i>^| and From Wick's theorem and the normal 

ordered form of Hamiltonian {Hn ^ H — (<i>„|i/|<i>„) = 

H — E^^^ ) , we get after the projection 

|i?jv + } + {-ff^2 } I = A£;f * 1 5i |$#0) 

{^z\hn + {h7si} + {h7s2}\^v) = AE:'\<i>z\s2m) 

In these equations, Ai?^** is the valence correlation en- 
ergy. It is defined as 



AEf, 



AE, 



N.corr 



^^JV-l,corr 



(12) 



FIG. 1: Diagrammatic representation of open shell cluster op- 
erators. The orbital lines with double arrows indicate valence 
and single up (down) arrow indicate particle (hole) states. 



where A£;^'™" and AE^-'^^'=°" are the total and core 
correlation energies respectively. The right members in 
Eq.( TopT l are what distinguishes the open shell coupled- 
cluster theory from that of the closed-shell. These are 
the equivalent of the folded diagrams in the many-body 
perturbation theory (MBPT) of open shell systems. 



A. Energy eigenvalue 

To obtain the energy eigenvalue E^ of the state \^y) 
project Eq. ^ on the state Then 



($„|i/(i + 5)|$„) = £;„, 



(13) 



here we have used (<i>i,|S'|<i>^,) — 0. Using the normal 
ordered Hamiltonian, defined earlier, Eq. ( 13 ) can be 
written as 



H 



N 



{l + S)\<^y)^Ey. 



From Wick's theorem 



Hn + {HnS} 



= AE., 



N,COTT 



(14) 



(15) 



The attachment energy is the difference in the exact en- 
ergy of the N- and {N — l)-electron state (closed-shell). 
In terms of correlation energies, attachment energy 



K = AE: 



-AE. 



N— 1 ,corr 



(16) 



where £„ is the single electron energy of the valence elec- 
tron. From the closed-shell coupled-cluster theory, the 
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FIG. 2: Diagrams which contribute to ABJ". The dashed 
line represent the residual Coulomb interaction. 

correlation energy A£;^"i'™"' have contribution fr om 
the closed diagrams. The right members in the amplitude 
equations Eq.( 10][TT | absorb this correlation energy as 
($:;;|{i?Ar5}|$t) is equivalent to AE^-'^'''°" {<^::.\S\<i>y) . 
Then the diagrams which contribute to AE'^** are the 
ones shown in Fig|2] 



B. Multiple valence shells 

It is relatively straight forward to calculate, from the 
single valence CCSD theory described, the ground state 
wave function and energy. Then the entire single particle 
basis space consist of one valence orbital, and the remain- 
ing are core (occupied) and virtual (unoccupied). How- 
ever, to calculate excitation energies, the excited atomic 
states and eigenvalues must be calculated. The trivial 
way is to solve the CCSD equations of each atomic states, 
ground and excited, separately. For example, to evaluate 
the 5c? excitation energy of Ba+ ion, the ground 

state |6s "^81/2) and the excited state \bd '^Dzj'il niust be 
calculated. Which translates to solving two sets of CCSD 
equations with ag^jBa^^) and |Ba^~'') as reference 



Here, iBa 



2+\ 



5^3/2 I 

is the closed shell Ba^+ reference 



states 
state. 

A better approach is to solve the ground and excited 
states CCSD equations in a single calculation. Then the 
theory is multi reference in nature and the cluster equa- 
tions of different states are coupled. In the present case, 
we choose the model space to consist of one state of spe- 
cific J and parity. Hence we do not have to invoke a full 
fledged multi reference coupled-cluster theory. 
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FIG. 3: Selected leading diagrams contributing to the hyper- 
fine structure constants in Eq.(21l. The dashed lines termi- 



nated with a circle represent hyperfine interaction. 



valence atoms, the magnetic dipole hyperfine structure 
constant 



Vj.(i« + l)(2j. + l) 



(19) 



Here, gi (/i = gjIfiN) is the gyromagnetic ratio and 
is the nuclear magneton. 



B. Hyperfine constants from coupled-cluster 



III. PROPERTIES CALCULATION 



A. Hyperfine Structure Constants 

The hyperfine interaction -ffhfs is the coupling of the 
nuclear electromagnetic moments to the electromagnetic 
field of the electrons. This causes splitting of the atomic 
levels and total angular momentum F is the conserved 
quantity. The atomic states are then \{IJ)FMf), here / 
and J are the nuclear spin and total electronic angular 
momentum respectively. The general form of the inter- 
action is 120 



(17) 



where t^{r) and T^' are irreducible tensor operators of 
rank k effective in the electron and nuclear spaces re- 
spectively. From the parity selection, only even and odd 
values of k are allowed for electric and magnetic inter- 
actions respectively. For the magnetic dipole interaction 
(/c = 1), the explicit form of the tensor operators are 



1 _ -^^/2[a■Cl(f)], 



and = iiq. (18) 



Here, Ci(r) is a rank one tensor operator in electron 
space and is a component of /x, the nuclear mag- 
netic moment operator. Then the nuclear moment is 
the expectation value of fj. in the stretched state = 
Parameters which represents the hyperfine 
splitting are the hyperfine structure constants. For one 



The measured value of an atomic property A for the 
atomic state is the expectation 



(20) 



In the present case, A is the hyperfine interaction i?hfs 
and in particular the magnetic dipole hyperfine interac- 
tion. From here on we use -ffhfsi however, the derivations 
and discussions are general, applicable to any dynami- 
cal variable. When coupled-cluster wave functions, from 
Eq.(|3]), are chosen as the correlated atomic states 

{^.,\Hms\^,) = {^,\Hms+2S^ Hms+S^ HmsS\^.), (21) 

where, -ffhfs = Hi^fse^ is the dressed operator. The 
factor of two in the second term on the right hand side 
accounts for HmsS as S'^ifhfs = H^isS. An expansion of 
Hhis ideal for an order wise calculation is 



n=l 



(22) 



The normalization factor, denominator in Eq.(20), in 
terms of coupled-cluster wave function is 



{^,\{1 + S^)e'^^e'^{l + S)\<^>,). (23) 



The dressed operator -ffhfs and operator in the 

normalization factor are non terminating series. In the 



next section we describe a method to calculate ffhfs to all 
order in T iteratively. To our knowledge, this is the first 
ever implementation of such a method within coupled- 
cluster theory. 



IV. PROPERTIES TO ALL ORDER 

For accurate properties calculations it is appropriate 
to include higher order terms in i?hfs- It is however non 
trivial to go beyond the second order, the number of di- 
agrams is large and a systematic evaluation is extremely 
tedious. On the other hand, diagrams can be grouped 
into different level of excitation (Zoe) and evaluate or- 
der wise iteratively. Here, loe is the number of core or 
valence electrons replaced with virtual electrons. For ex- 
ample, the diagrams in Fig|4] have loe one. In each of 
these diagrams, one core electron is replaced by a virtual 
electron. 

To calculate the diagrams of loe one to all order, con- 
sider the loe one diagrams arising from i/hfse"^. That 
is 



'hfs 



HmsT 



1 



(24) 



Where the subscript denotes the loe of the contributing 
terms. It is equivalent to a one-particle interaction and 
considered as effective properties operator which incor- 
porates electron correlations. In the next iteration 



1 



(25) 



where i = 1 , 2 in CCSD and the superscript conn imply 
only the connected diagrams contribute. From the defi- 
nition of the cluster operators, Ti and Tj have loe i and 
—i respectively. The above equation is equivalent to the 



expression in Eq.(24l sandwiched between cluster oper- 
ators of equal but opposite loe. So the net loe remains 
unchanged. In general, we can then write 

(26) 

This is an iterative equation and it is possible to evalu- 
ate it order by order to convergence. The sum of all the 
contributions is equivalent to calculating the effective op- 
erator 



■Hi 



(27) 



This contribute to the hyperfine structure as 5*2 Tii. At 
the lowest level there are diagrams and correspond to 
Fig|^-k. In a similar same way, the effective properties 
of higher loe are calculated. 




FIG. 4: Diagrammatic representation of the iterative equation 
to calculate the loe one effective hyperfine operator Hmsi^ ■ 
The iteration is implemented with the and T2. 



For further study, we resort to diagrammatic analy- 
sis. Consider diagrams arising from {Hms&^)i, there are 
six diagrams in total. These are shown on the first row 
at the right hand side of Fig|4j These define the ini- 
tial choice of the effective diagram. For the next and 
higher iterations, consider the contractions with and 
T2- The contribution from the t\ and Ti is neglected 
as these cluster amplitudes, on an average, are several 
orders of magnitude smaller than T2. Then the iteration 
is equivalent to the diagrammatic equation in Fig|4] and 
it is mathematically 



(28) 



Where Ti!l is (iJhfse"^)i, the effective operator prior to the 
iteration. Since only the unique diagrams are considered, 
there are no multiplying factors. The algebraic relation 
in Eq.(28l is also without multiplying factors as the se- 



quence of the contraction is uniquely defined. Which is 
not the case in the expansion of e^^^iJhfse"^. 



V. DESCRIPTION OF NUMERICAL METHODS 

The calculations presented in the paper involve vari- 
ous numerical techniques and methods. Some are fairly 
straight forward, oft used in atomic theory calculations. 
Others are not, specialized and application specific. For 
easy reference in future works, we provide an outline of 
the numerical methods used. This is appropriate as we 
recommend, based on the current work, an approxima- 
tion of the properties operator in coupled-cluster theory. 



A. Atomic Hamiltonian and single particle states 

In the results presented in this paper the Dirac- 
Coulomb Hamiltonian is chosen H^*^ for the calculations. 
It incorporates relativity at the single particle level accu- 
rately. And, as the name indicates, the Coulomb inter- 
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actions between the electrons. For an N electron atom 



B. Basis set and cluster amplitudes 



H 



N 



t.-p. + ip-iy-VNir,)] + 



E 

l<j 



(29) 

where p is the linear momentum, and and (3 are the 
Dirac matrices. For the nuclear potential VN^r), we con- 
sider the finite size Fermi density distribution 



Pnuc('') = 



Po 

1 + e^^^'^^"- ' 



(30) 



here, a — t41n3. The parameter c is the half-charge 
radius, that is /9nuc(c) = /5o/2 and t is the skin thickness. 
At the single particle level, the spin orbitals are of the 
form 

^ Pnnir)xnm{r/r) \ 

Pn.m{r) ~ ^ tQ,,,{r)x-.Ar/r) J ' ^^^^ 

where PnK,{f) and Qn^if) are the large and small com- 
ponent radial wave functions, k is the relativistic total 
angular momentum quantum number and x„m(r/r) are 
the spin or spherical harmonics. One representation of 
the radial components is to define these as linear com- 
bination of Gaussian like functions and are referred to 
as Gaussian type orbitals (GTOs). Then, the large and 
small components [211 ES] are 



(32) 



p 

Q„K(r) = 5]cfp5fp(0- 



The index p varies over the number of the basis functions. 
For large component we choose 

g^p{r)^C^.,r"^e-'^^^\ (33) 

here is an integer. Similarly, the small component are 
derived from the large components using kinetic balance 
condition. The exponents in the above expression follow 
the general relation 



(34) 



The parameters ao and /3 are optimized for each of the 
ions to provide good description of the properties. In 
our case the optimization is to reproduce the numerical 
result of the total and orbital energies. The optimized 
parameters used in the calculations are listed in Table |T] 
From Eq.(31 1 the reduced matrix element of the mag- 
netic hyperfine operator between two spin orbitals , v' 
and V, is 

dv 

A detailed derivation is given in Ref. [23]. 



For all the alkaline Earth metal ions considered, Mg+, 
Ca^, Sr+ and Ba^, we use V^~'^ orbitals for the calcula- 
tions. This is equivalent to calculating the spin orbitals 
from the single particle eigenvalue equations of the dou- 
bly ionized alkaline Earth metal atoms, namely Mg^+, 
Ca^+, Sr^+ and Ba^+. Then the single particle basis 
sets have few bound states and rest are continuum. We 
optimize the basis such that: single particle energies of 
the core and valence orbitals are in good agreement with 
the numerical results. For this we use GRASP92 [53] to 
generate the numerical results. 

As mentioned in earlier sections, we compute the 
closed- shell cluster amplitudes T first. These are used 
to generate the open shell cluster amplitudes S. The 
coupled nonlinear and linear equations are solved itera- 
tively. We employ direct inversion in the iterated sub- 
space (DIIS) ^25j for convergence acceleration. 



VI. RESULTS AND DISCUSSIONS 

A. Ionization potential and excitation energies 

To determine the quality of the basis set and parame- 
ters, we compute the attachment energies of the ground 
state (51/2) and the first excited P1/2, P3/2, D3/2 and 
states are calculated. Then the ionization potential 
(IP), the energy required to remove the valence electron, 
is the negative of the attachement energy —E'^^^. To 
calculate the excitation energy (EE) of the state \'9v), 
consider Eg^^ and E^^^ as the attachment energies of 
the ground state and excited state. Then difference 
jjatt _ ^att -g ^Yic EE, it can as well be defined in terms 
of IPs. 

For further analysis on the correlation effects incor- 
porated with CCSD, we first compute IP with relativis- 
tic many-body perturbation theory (MBPT). The MBPT 
calculations are similar to our previous work [26] for sec- 
ond order correlation energy of closed-shell atoms, in par- 
ticular noble gas atoms. The MBPT diagrams of IP are 
similar to the first four attachment diagrams in Fig|2]but 
have residual interaction in stead of S and T operators. 
Where, the first two diagrams Figj2]^a-b), direct and ex- 
change, have the valence replaced by a virtual state and 
encapsulates core- valence correlation. The remaining di- 
agrams Fig|2|c-d) represent core-core correlation as these 
involve double replacement of core electrons. The other 
two diagrams with Ti do not contribute as single replace- 
ments with residual Coulomb interaction are zero. 

The results of the MBPT calculations are listed in Ta- 
ble|lll The 3d ^^^3/2,5/2 and 3p ^Pi/2.,3/2 of Mg+, evalu- 
atedirom the MBPT IPs, are marginally lower than the 
experimental data but are very close. From Ca"*", there 
is a change in the pattern of the EEs. The MBPT results 
of ^£'3/2,5/2 EEs are lower than the experimental data, 
whereas the ^^1/2. 3/2 EEs are higher. The same pattern 



6 



occurs in Sr+ and Ba+. Similar pattern is observed in 
the results of previous calculations [28]. The differences 
between the results in Ref.[55] and ours are minor and 
random in nature. These deviations can be attributed to 
the nature and completeness of the basis sets chosen in 
the two calculations. 

The CCSD results of the EE are also listed in Table|!l] 
these are closer to the experimental data than the MPBT 
results. This is not surprising as CCSD encapsulates 
electron correlations more accurately. The trend of the 
CCSD results separates into two: Mg"*" and other ions. 
The additional electron correlation increases the IPs of 
Mg+ , whereas there is a decrease in the IPs of Ca+ , Sr+ 
and Ba+. However, the states change differently such 
that the EE improves. These results vouch for the relia- 
bility of the basis set for properties calculations. 



B. Magnetic dipole hyperfine constants 

To compute the hyperfine constants from the CCSD 
wave functions, we use Eq.(21l. The results are listed 



in Table |III| for comparison the results of other theoret- 
ical calculations and experimental data are also given. 
As defined in Eq.(21l, the coupled-cluster expression of 



the hyperfine structure constants is separated into three 
groups. The dominant contribution from the first term 
-^hfsi up to first order in and T, is 



Hhfs ~ Hi^fs + 2-ffhfsT'i 



2Tn 



(36) 

Here, the first term is the Dirac-Fock (DF), which has the 
largest contribution. The factor two in the second and 
fourth terms accounts for the complex conjugate terms. 
The third term, second order in Ti, has one diagram and 
negligibly small contribution. The diagrams arising from 
the last term are topologically similar to the attachment 
diagrams (c-d) in Figl2| However, with the in stead 
of residual Coulomb interaction and i?hfs inserted on one 
of the orbital lines. There are ten diagrams and contri- 
bution from these are labelled as i/hfs — DF. The last 



two terms in Eq.(21l are approximated as 



S^HhfsS 



2S^ (i/Mse^)i, 
sIHms{Si+2S2) 



(37) 

sIHmsS2- (38) 



Like in ifhfs , the factor of two is to account for the com- 
plex conjugate terms. The expression of (7?hfse^)i is as 
given in Eq.(24l. The S'|i?hfsS'2 term have contributions 
from the diagrams (b-gHn Figj3] These are topologically 
similar to (a-b) in Figji] But, like in rj'iJhfsTb, Sl in- 
stead of residual Coulomb interaction and i?hfs inserted 
to one orbital line. Diagrams arising from the remaining 
terms are also given in Fig|3] Based on this grouping, 
the contributions are listed in Table |IV] In the following 
we present a detailed comparison of our magnetic hyper- 
fine constants results with the earlier ones. As discussed 



later, some of our results are the best match with ex- 
perimental data. This is a thorough test for the starting 
point of our iterative procedure and the expression for 
properties calculation we recommend. 



1. Mg+ 

The experimental data is available only for the ground 
state 3s [S]- However, theoretical results are avail- 

able for the low lying states 3s '^Si/2, 3p '^Pi/2, 3p ^P^/2, 
3d ^D^/2 E^ud 2>d ^D^i2- In the previous works, the calcu- 
lations used relativistic many-body perturbation theory 
[29 , 30 and linearized CCSD ^31| using numerical and B- 
splines basis sets respectively. These report the DF con- 
tribution for 3s 2S'i/2 as -466.4 P] and -463 [23]. The 
later is in excellent agreement with our result —463.29. 
The other dominant terms are S'^i/hfs and iJhts — DF, 
contribution from these are —107.32 and —16.13 respec- 
tively. Total value of these three terms is —586.76, 98% 
of the experimental value. Our total value —596.78, after 
including S^HmsS, is 0.08% lower than the experimental 
value and is the best theoretical result. 

For the 3p ^Pi/2 and 3p '^Pz/2 states, our DF val- 
ues —76.98 and —15.24 are in very good agreement with 
the values —77 and —15.2 given in Ref. [29 . The 
two states have 20.6% and 16.3% contribution to the 
total value from S^H^fs- This difference shows varia- 
tions in the nature of correlation effects, predominantly 
core-polarization. Our total values for the two states are 
— 103.0 and —19.55, these are in very good agreement 
with the previous results. 

For the id'^D^/2 and id'^D^i2 states, our DF values 
are —1.26 and —0.54 respectively. Whereas, the values 
in a previous work [55] are —1.61 and —0.54. The results 
of 3d^D^/2 match perfectly but there is a significant dif- 



ference in the results of 3d D 



3/2- 



Our result of -1.26 



is 28% less in magnitude. The correlation effects, core- 
polarization in particular, are markedly different from 
the other states. Contribution from 5^i/hfs to 3d'^D^/2 
is 0.19 is 15% in magnitude of Dirac-Fock and opposite 
in sign. It is even more dramatic for 3d^Dc^/2-, it is 0.65, 
which larger than Dirac-Fock in magnitude and opposite 
in sign. 

Considering that the calculations in Ref. [25] incorpo- 
rates core-polarization to all orders , we can extract the 
pair correlation effects. For the 3s ^6*1/2 state, the core- 
polarization contributes —91. Subtracting this from our 

Hhfs result, the pair correlation contribution to this 
term is —16.32. Adding the other terms as well, the to- 
tal contribution from pair correlation is —43.82. Which 
is less than the core-polarization but not negligible. For 
the other states the core-polarization contributions are 
-18, -3.7 and 0.71 for 3p '^Pi/2, 3p ^Pg/a and 3d ^Dg/a 
respectively. The corresponding pair correlation contri- 
bution are —7.86, —0.62 and —0.04. The pair correlation 
is negligible in last two states and we have not estimated 
for 3d ^D^/2 as there is a large difference between our 



7 



DF value and Ref. ^ . 



2. Ca+ 

This is the most well studied, experimentally and the- 
oretically, singly ionized alkaline Earth ion. There is a 
large variation in the experimental results of 4s ^5*1/2 and 
Ap and less in the results of 4p ^^3/2, Sd ^D^f2 

and 3d states. On the other hand the theoreti- 

cal results exhibit significant variations for all the states 
except 4p ^^3/2. The DF values of 4s reported in 

previous works are —589 [29] and —588.933 [32], these are 
calculated with numerical and B-spline basis sets respec- 
tively. Our value —589.09 is in very good agreement with 
these results. The core- valence correlation from S^H\ifs 
accounts for 22% of the total value. This is much larger 
than in Mg+ (17%). On the other hand core-core corre- 
lation, contribution from i?hfs — DF, is smaller. Our total 
value —808.12 is marginally higher than the experimental 
values but lies between the other theoretical results. 

In previous studies DF values of the Ap ^Pi/2 are —102 
[29] and -101.492 |32]. Similarly, for 4p ^P3/2 the val- 
ues are -19.2 [55] and -19.646 [32]. These are in very 
good agreement with our results —101.47 and 19.65. Like 
in 4s ^5i/2 there is an increase, compared to Mg+, in 
S^Hhfs contribution. It accounts for 26% and 30% of 
the total value for the two states. Our total values 
of 4p is lower than the other theoretical results. 

Whereas ip ^^-3/2 exhibits opposite trend. 

For 3d ^0^/2, the DF values in the previous studies 
are -33 [IS], -33.206 [32 and -39.12 [33.. The first 
two compares well with our value —33.55. Similarly, our 
3d ^1)5/2 DF value —14.29 is in good agreement with the 
previous results -14 p9] and -14.144 [32]. There is a 
change in the nature of 5^ifhfs contribution to 3d '^0^/2- 
Unlike in Mg+, it is in phase with DF and similar trend 
is observed in Sr^" and Ba"*" as well. The contribution 
from {S^Hhfs + c.c.) to 3d ^1^5/2 is the only one which 
is less in magnitude than the DF value. In all the other 
ions ( Mg+, Sr+ and Ba+) DF values are less in mag- 
nitude. The impact of core-core correlation is not large 
but not negligible. Our total value for 3d ^I?3/2 is lower 
than all the theoretical and experimental values. How- 
ever, our result for 3d ^D^j2 matches very well with the 
experimental data. 

Taking the core-polarization results from Ref. and 
following the procedure in Mg+ we estimate the pair cor- 
relation effects. We get the pair correlation contributions 
as -108.61, -19.37, -4.25, -11.96 and -7.99 for the 
4s 2^1/2 and Ap ^Py2, 4p ^Pg/^ 3d 2^^3/2 and 3d ^0^/2 
respectively. Except for 4p 2P3/2 and 3d 2-D3/2, these are 
in very good agreement with the pair correlation listed 
in Ref. [32] • Not surprisingly, our results for these two 
states deviate from the other theoretical and experimen- 
tal data. 



3. Sr+ 

Experimental data is limited to 5s '^81^2, Ap ^Pi/2 
and 4d ^D^/2. However, several theoretical investi- 
gations have examined the hyperfine structure of Sr^. 
The 5s "^81/2 DF value earlier works are —735 [34' and 
-736.547 [55]. Our value -738.204 is higher than both 
of the values. There is a large contribution from S'^i/hts- 
It is 22% of the total value and same as 4s "^81/2 of Ca+. 
The core-core correlation is less significant. Our total 
result is lower than the experimental data and other the- 
oretical results. 

The DF values of bp '^P\/2 from previous works are 
-122 [34j and -121.576 [32]. And values for bp 2P3/2 
are -21.4 [51] and -21.331 [52]. These are less than 
our values —122.363 and —21.501. The core-core corre- 
lation effects is negligibly small, 0.3% of the total value. 
Compared to Ca+ {Ap ^Pj), there is an enhanced role 
of S^Hhfs in bp 2^3/2- It accounts for 33% of the total 
value. Our total value for bp '^Pi/2 is less than the pre- 
vious theoretical results. But, the value of bp ^Pi/2 is in 
excellent agreement with the experimental data. 

The DF values of 4d ^0^/2 from previous works are 
-31.2 34J, -31.126 [52] and -34.23 [55]. And values 
for 4d '^D^/2 are -13.0 [5T, -12.977 [5T and -14.27 
[33] . These compare well with our values —31.368 and 
— 13.080. There is a marked change in the role of S^Hhfs 
for the 4d '^D^/2 state. It has larger magnitude ( 135% ) 
than the DF value. Our total value for 4d ^Dy,/2 is lower 
than the other theoretical values. However, 4d ^1)3/2 is 
in excellent agreement with the experimental data. 

There is noticeable difference in the estimates of the 
core-polarization effects in the earlier works [521153]. For 
example, the core-polarization contribution to 4d "^Dzji 
is estimated as —6.3 in Ref. ,34J, whereas it is —2.413 in 
Ref. [52]. For consistency of analysis, with the choice in 
Ca"*", we adopt the core-polarization results of Ref. [52] 
and estimate pair correlation effects in our results. These 
are -127.795, -24.241, -3.938, -11.235 and -6.594. Af- 
ter accounting for the difference in the Dirac-Fock results, 
the results for 5s '^8^12 and 4d ^D^/2 are in very good 
agreement with Ref. [32] . 

4. Ba+ 

It is a candidate system, as mentioned earlier, for 
a novel technique to measure parity nonconservation 
(PNC) experiment [TS]. In this context, theoretical study 
of Ba+ hyperfine constants is very important. It is a good 
proxy for the PNC in ions or atoms arising from neutral 
weak currents. Except for 6p '^Pi/21 there are experimen- 
tal data for the low-lying states and theoretical results 
are available for 6s "^81/2, 6p ^Pi/2, 6p 2/3/2, 5d ^Dg/j 
and 5d '^D^i2- The DF value of 6s "^8x12 in previous 
calculations are 2929.41 f5F and 3055 [51]. Our result 
is 3003.105 and significantly different from both of the 
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values. The contribution from the core-core correlation 
Hilts — DF is of opposite phase to the DF contribution. 
This is in contrast to the states we have discussed so 
far. The total value is in very good agreement with ex- 
perimental data but significantly different from the other 
theoretical results. It has 23% contribution from Hms- 

The DF value of 6p and 6p in the previous 

work are 492.74 ^ and 71.84 |3l]. These are different 
from our values of 504.196 and 73.674. The core-core cor- 
relation i/hfs — DF, as in 6s ^51/2, is of opposite phase for 
6p The total results of the two states are 705.039 

and 130.191. The first is lower than the previous theoret- 
ical result. And the later is in very good agreement with 
the theoretical result but lower than the experimental 
data. 

The DF values of 5d ^1)3/2 in the previous studies are 
128.27 [35] and 139.23 [33] • And for hd ^1)5/2 the values 
are 53.213 [37] and 55.82 [33]. Our results are 129.875 
and 52.085, these are closer to Ref.[35] and Ref.[37] re- 
spectively. The Hms contribution to ^i?5/2 is large, 
141% of the DF value and opposite in phase. Our total 
total values 185.013 and —12.592 are close to experimen- 
tal data. 

For Ba+, except for the ground state there are no sys- 
tematic studies of core-polarization effects. The previous 
work of Ref. [35] uses methods, GTO basis and relativis- 
tic coupled-cluster, similar to what we have used in the 
present paper. Comparing the two, there is a good corre- 
lation between different coupled-cluster terms for all the 
states except 6s ^5*1/2. 



All order calculations 



tized representation of H.. That is 



ijkl 



Tna\.ai ajQi + ■ ■ ■ 



(39) 



The Eq.(28l then assumes the form 



qr*qrrpr , Pf*qrrqr r,*qrlqp , .r.-. 

H ''be ^ca + '^bc ''ah + ''be ''be ' \^^) 

where hij is the matrix element (i|ft.hts|j) and Vf^ = t'^!- — 
t'jl is the antysymmetrised cluster ampHtude. This is the 
equation we solve iteratively till convergence. After each 
iteration, we evaluate the contribution from the effective 
operator to the hyperfine constant 5|7ii. The results of 
our calculations are given in Table fV] For most of the 
cases, the results converges to KHz accuracy after two 
iterations. 

In terms of absolute changes, the largest is observed 
in 6s ^5i/2 of Ba+. For this state the zeroth iteration, 
arising from Tii, as given in Tablejv] is 469.636. It con- 
verges to 467.450 at the second iteration and change is 
-2.186. Which is 0.5% of the zeroth iteration and 0.05% 
of the total value. Whereas in terms of fractional change, 
the largest is 5d Upon convergence the change is 

-0.702, which is 1.9% of zeroth iteration. However, this 
is 5.5% of the total value. This arises from the large 
cancellation between the DF and S^H^ifg. Here to obtain 
correct result the iterated calculation should be applied 
to the other terms as well. Not very surprisingly, the 
changes in Mg+, Ca+ and Sr+ which have lower Z are 
negligibly small. 

Considering that iteration is implemented for the loe 
which contributes maximally. Contributions from the 
other loe is expected to much smaller. 



In the previous section we analyzed and compared our 
results with the earlier ones in some detail. Majority of 
our results are in very good agreement with the exper- 
imental data, some are perhaps the best match. The 
earlier works chosen for comparison are based on di- 
verse types of orbitals: numerical, finite discrete spec- 
trum, B-spline and GTO. These are an accurate repre- 
sentation of the tried and tested types of single orbital in 
atomic calculations. Similarly, there is a variation in the 
many-body methods: MBPT, MCDF-EOL and coupled- 
cluster. This is a large data set for comparison. More 
importantly, among the ions considered there are large 
changes in the role of electron correlations. This choice is 
essential to deliberate on the consequence of higher order 
terms and avoid erroneous inference from an incomplete 
sample. This sets the stage for a systematic appraisal of 
the higher order terms in the properties calculations. 

As discussed in Section |IV] we implement the iterative 
method to calculate the hyperfine constant to all orders 
for the loe one. To frame the iterative equation in terms 
of components, define r as c numbers in the second quan- 



VII. CONCLUSIONS 

We have calculated, as well as surveyed and compared 
the magnetic hyperfine structure constants of low lying 
states of Mg+, Ca"*", Sr+ and Ba"*" available in the litera- 
ture. For the states 3s "^81/2 (^^Mg+), 3d ^D^/^ C'^Ca+), 
4d (^^Sr+) and 6s 251/2 ("^Ba+), our results pro- 

vides the best match with the experimental data. Further 
more, results for most of the other states are in very good 
agreement with the available experimental data. 

The chosen systems have hyperfine constants with 
varying dependence on electron correlations. It is a suit- 
able data set to examine the impact of higher order terms 
in properties calculations with relativistic coupled-cluster 
theory. This is of paramount importance for high pre- 
cision properties calculations with relativistic coupled- 
cluster. Our study establish without any ambiguity, the 
higher order terms are not important when the leading 
terms DF and S^Hms contribute coherently. However, 
when large cancellation occurs like in 2-D5/2 state of al- 
kaline Earth ions, a consistent calculation of the different 
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terms to equal orders is a must. Except for such cases, 
based on the present study, we recommend 

(*„|Hhfe|*^) = ($,|Hhfe + 2i/hfeTi+T//fhfs(ri + 2T2) 

+SlHhts {Si + 2S2) + SlHhi,S2\^^}i41) 

to calculate hyperfine and similar properties for single 
valence systems. It is sufficient to include terms up to 
quadratic in T for properties c;alc;ulations. Higher order 
terms, all together, have less than 0.1% of the total value 



and can be neglected. 
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TABLE I: Optimized parameters a and /3 of the GTO basis used in the calculations. 



Symmetry 




25Mg+ 






«Ca+ 






«^Sr+ 






"^Ba+ 






a 


/3 


Basis 


a 


/3 


Basis 


Q 




Basis 


a 





Basis 








function 






function 






function 






function 


s 


0.0083 


2.8900 


28 


0.0063 


2.8800 


29 


0.0083 


2.9800 


30 


0.0063 


2.9800 


31 


V 


0.0072 


2.9650 


25 


0.0072 


2.9650 


26 


0.0072 


2.9650 


27 


0.0072 


2.9590 


28 


d 


0.0070 


2.7200 


22 


0.0070 


2.7000 


24 


0.0070 


2.8000 


25 


0.0070 


2.4500 


26 



TABLE IL Ionization potential and excitation energies. For comparison other results and experimental values are also listed. 
All values are in atomic units. 



Ion 


state 


MBPT 


Coupled-cluster 


Other works 


Exp results Rcf 27 . 






IP 


EE 


IP 


EE 


IP 


EE 


EE 


25Mg+ 


3Sl/2 


-0.55156 


0.0 


-0.55203 


0.0 


-0.55252 


0.0 


0.0 




3d3/2 


-0.22652 


0.32504 


-0.22666 


0.32537 


-0.22677 


0.32575 " 


0.32573 




3^5/2 


-0.22652 


0.32504 


-0.22668 


0.32535 


-0.22677 


0.32575 " 


0.32574 




3pi/2 


-0.38922 


0.16234 


-0.38950 


0.16253 


-0.39003 


0.16249 " 


0.16252 




3P3/2 


-0.38878 


0.16278 


-0.38917 


0.16286 


-0.38961 


0.16291 " 


0.16294 


«Ca+ 


4Si/2 


-0.43784 


0.0 


-0.43671 


0.0 


-0.43836 


0.0 


0.0 




3(^3/2 


-0.37797 


0.05987 


-0.37601 


0.06070 


-0.37768 


0.06068 ' 


0.06220 




3d5/2 


-0.37762 


0.06022 


-0.37578 


0.06093 


-0.37731 


0.06205 ' 


0.06247 




4pi/2 


-0.32180 


0.11604 


-0.32128 


0.11543 


-0.32217 


0.11619 ' 


0.11478 




4P3/2 


-0.32075 


0.11709 


-0.32119 


0.11552 


-0.32111 


0.11725 


0.11580 


«^Sr+ 


5Sl/2 


-0.40788 


0.0 


-0.40573 


0.0 


-0.40839 


0.0 


0.0 




4^3/2 


-0.34236 


0.06552 


-0.33926 


0.06647 


-0.34279 


0.06560 ' 


0.06632 




4^5/2 


-0.34091 


0.06697 


-0.33827 


0.06746 


-0.34132 


0.06707 '• 


0.06760 




5pi/2 


-0.29793 


0.10995 


-0.29696 


0.10877 


-0.29838 


0.11001 ' 


0.10805 




5p3/2 


-0.29421 


0.11367 


-0.29425 


0.11148 


-0.29463 


0.11376 ' 


0.11171 




6Sl/2 


-0.37297 


0.0 


-0.36862 


0.0 


-0.37308 


0.0 


0.0 




5^3/2 


-0.35296 


0.02001 


-0.34758 


0.02104 


-0.35172 


0.02136 ' 


0.02221 




5^5/2 


-0.34872 


0.02425 


-0.34386 


0.02476 


-0.34748 


0.02560 ' 


0.02586 




6pi/2 


-0.27685 


0.09612 


-0.27483 


0.09379 


-0.27532 


0.09776 


0.09232 




6p3/2 


-0.26882 


0.10415 


-0.26821 


0.10041 


-0.26946 


0.10362 ' 


0.10002 



"ReferencelSl]. 
''Reference[28]. 
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TABLE III: Magnetic dipole hyperfine structure constants (in MHz) for ''Mg+, ■*''Ca+, **^Sr+, and "^Ba+ ions 



Ion 


state 


This work 


Other works 


Experiment 


^5]V[g+ 


3Si/2 


-596.785 


-597.6', -554", -(602 ±8)' 


-596.254™ 




3pi/2 


-102.997 


-103.4\-100" 


- 




3p3/2 


-19.546 


-19. 29', -19" 


- 




3(^3/2 


-1.083 


-1.140', -1.25" 






3^^5/2 


0.118 


0.1196', 0.107'',0.17" 




'*^Ca+ 


4Sl/2 


-808.126 


-805.35\ -819», -794.7'',-806.4(2.5)'' 


-797.5(2.4)" -805(2)'' 




4pi/2 


-142.782 


-143.07*,-148», -144. 8'', -143", -145.4(4)" 


-158(3.3)", -145.5(1.0)'', -142(8)", -145.4(0.1)^ 




4P3/2 


-32.185 


-30.50*, -30.9», -29.3'',-30", -30.4(4)" 


-29.7(1.6)", -31.9(0.2)'', -31.0(0.2)-^ 




3^3/2 


-45.294 


-47.82', -523,-49.4'', -47.3(3)" 


-48.3(1.6)", -47.3(0.2)^ 




"-"^5/2 


—4.008 


—3 351" —3 55* —5 2^ —4 2'' —3 6f3~)" 


—3 8(0 gV 3 893^21" 




5Si/2 


-990.638 


-10003.18\ -1000* 


-1000.5(1.0)' 




5pi/2 


-169.988 


-178.40^,-177*, -175" 


- 




5P3/2 


-36.225 


-35.11^,-35.3* -30" 


-36.0' 




4^3/2 


-44.320 


-47.36', -46.7* 






4d5/2 


2.168 


2.156 , 2.51 , 1.1 


2.1 


"7Ba+ 


6Si/2 


4021.721 


4072. 83" 


4018.2' 




6pi/2 


705.039 


736.98" 






6P3/2 


130.191 


130. 94", 126" 


126.9°, 112.77'' 




5^3/2 


185.013 


192. 99", 188. 76" 215° 


189.730°, 170.88'' 




5^5/2 


-12.593 


9.39",-11.717", -18° 


-12.028° 



"Reference 37j 
'Referencegg 
"ReferencefSH] 
''Reference 39' 
'^Reference 40 j 
.^Referencel^l 



5 Reference 
''Reference 
'Reference 
^Reference 



*^ Reference l3l 
& 

M 
H 

25, 

46 



'Reference 
'"Reference 
"Reference 
"Reference 
^Reference 
5 Reference 
""Reference f47 
^Reference 129 
'Reference :30 
"Reference Ho 
■"ReferenceES' 
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TABLE IV: Contributions from different terms in the coupled-cluster, magnetic dipole hyperflne constant, properties expression. 
The values listed are in MHz. 



ion 


state 








Coupled-cluster terms 












DF 


- 

Hmb-'DF 




cf TT C 


Ot TT CI 

y3{/lhfs>3l 


TT C 




Norm 










+C.C 


+C.C. 


























+C.C. 






3Sl/2 


-463.297 


-16.136 


-107.325 


-1.532 


-0.396 


-5.560 


-3.878 


1.002 




3pi/2 


-76.984 


-2.754 


-21.254 


-0.326 


-0.089 


-0.989 


-0.720 


1.001 




3p3/2 


-15.242 


-0.695 


-3.184 


0.005 


-0.018 


-0.277 


-0.160 


1.001 




3c?3/2 


-1.259 


-0.007 


0.186 


0.004 


-0.001 


-0.008 


-0.000 


1.001 




3(^5/2 


-0.540 


-0.003 


0.648 


0.017 


-0.000 


-0.004 


-0.000 


1.001 


''^Ca+ 


4si/2 


-589.087 


-12.696 


-180.217 


-4.313 


-1.802 


-10.717 


-16.771 


1.009 




4pi/2 


-101.473 


-0.497 


-37.514 


-0.978 


-0.446 


-1.031 


-1.638 


1.006 




4P3/2 


-19.648 


-0.321 


-9.691 


-0.222 


-0.094 


-1.004 


-1.426 


1.007 




3(^3/2 


-33.554 


-2.553 


-7.701 


-0.149 


-0.260 


-2.153 


0.258 


1.018 




3(^5/2 


-14.294 


-1.247 


13.430 


0.481 


-0.111 


-2.449 


0.110 


1.018 


*''Sr+ 


5Sl/2 


-738.204 


-3.667 


-218.305 


-5.258 


-3.046 


-15.027 


-18.379 


1.011 




5pi/2 


-122.363 


-0.675 


-44.231 


-1.120 


-0.678 


-1.446 


-0.637 


1.007 




5P3/2 


-21.501 


-0.398 


-12.011 


-0.339 


-0.126 


-1.043 


-1.099 


1.008 




4(^3/2 


-31.368 


-3.084 


-8.431 


-0.271 


-0.139 


-1.979 


0.255 


1.016 




44/2 


-13.080 


-1.626 


17.644 


0.470 


-0.058 


-1.991 


0.843 


1.016 


"7Ba+ 


6Sl/2 


3003.105 


-39.093 


939.272 


23.989 


17.598 


66.108 


68.032 


1.014 




6pi/2 


504.196 


-5.948 


196.073 


5.397 


4.072 


6.064 


1.982 


1.010 




6p3/2 


73.674 


0.665 


45.835 


1.555 


0.619 


4.480 


4.797 


1.011 




54/2 


129.875 


12.565 


37.918 


1.148 


0.462 


9.495 


-2.329 


1.022 




54/2 


52.085 


7.240 


-73.611 


-1.520 


0.191 


9.554 


-6.803 


1.022 
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TABLE V: Magnetic dipole hyperfine structure constant, contributions from higher-order terms in the all order scheme Eq. ( 22 1 



I n 

on 








J -nhfs 










ILCi — U 


iter 1 


ifor — 9 

iLei — Z 


ix6r 


Oonverged 






(7?hfse^)i 


rj'(i?hfse^)ir2 


TJ^//hfte^)ir2^ 


Tf{HM.e^)iTi 


value 


25Mg+ 


3Si/2 


-53.663 


-53.502 


-53.503 


-53.503 


-53.503 




3pi/2 


-10.627 


-10.563 


-10.564 


-10.564 


-10.564 




3p3/2 


-1.592 


-1.577 


-1.577 


-1.577 


-1.577 




3(^3/2 


0.093 


0.091 


0.091 


0.091 


0.091 




3t?5 /2 


0.324 


0.321 


0.321 


0.321 


0.321 


^^Ca+ 


4Sl/2 


-90.109 


-89.776 


-89.778 


-89.778 


-89.778 




4j5l/2 


-18.757 


-18.570 


-18.574 


-18.574 


-18.574 




4p3/2 


-4.845 


-4.792 


-4.793 


-4.793 


-4.793 




3^3/2 


-3.851 


-3.887 


-3.885 


-3.885 


-3.885 




3^5/2 


6.715 


6.638 


6.639 


6.639 


6.639 


^■^Sr+ 


5Si/2 


-109.153 


-108.716 


-108.720 


-108.720 


-108.720 




5j5l/2 


-22.116 


-21.908 


-21.912 


-21.912 


-21.912 




5P3/2 


-6.006 


-5.943 


-5.944 


-5.944 


-5.944 




4^3/2 


-4.216 


-4.267 


-4.265 


-4.265 


-4.265 




4(^5/2 


8.822 


8.687 


8.689 


8.689 


8.689 


"7Ba+ 


6Si/2 


469.636 


467.423 


467.450 


467.449 


467.449 




6pi/2 


98.036 


97.052 


97.075 


97.074 


97.074 




6P3/2 


22.917 


22.655 


22.660 


22.660 


22.660 




5^3/2 


18.959 


19.161 


19.150 


19.150 


19.150 




5(^5/2 


-36.806 


-36.092 


-36.104 


-36.104 


-36.104 



